Projecting the development of antimicrobial resistance in Neisseria gonorrhoeae from antimicrobial surveillance data: a mathematical modelling study.

BACKGROUND: The World Health Organization recommends changing the first-line antimicrobial treatment for gonorrhoea when ≥ 5% of Neisseria gonorrhoeae cases fail treatment or are resistant. Susceptibility to ceftriaxone, the last remaining treatment option has been decreasing in many countries. We used antimicrobial resistance surveillance data and developed mathematical models to project the time to reach the 5% threshold for resistance to first-line antimicrobials used for N. gonorrhoeae. METHODS: We used data from the Gonococcal Resistance to Antimicrobials Surveillance Programme (GRASP) in England and Wales from 2000-2018 about minimum inhibitory concentrations (MIC) for ciprofloxacin, azithromycin, cefixime and ceftriaxone and antimicrobial treatment in two groups, heterosexual men and women (HMW) and men who have sex with men (MSM). We developed two susceptible-infected-susceptible models to fit these data and produce projections of the proportion of resistance until 2030. The single-step model represents the situation in which a single mutation results in antimicrobial resistance. In the multi-step model, the sequential accumulation of resistance mutations is reflected by changes in the MIC distribution. RESULTS: The single-step model described resistance to ciprofloxacin well. Both single-step and multi-step models could describe azithromycin and cefixime resistance, with projected resistance levels higher with the multi-step than the single step model. For ceftriaxone, with very few observed cases of full resistance, the multi-step model was needed to describe long-term dynamics of resistance. Extrapolating from the observed upward drift in MIC values, the multi-step model projected ≥ 5% resistance to ceftriaxone could be reached by 2030, based on treatment pressure alone. Ceftriaxone resistance was projected to rise to 13.2% (95% credible interval [CrI]: 0.7-44.8%) among HMW and 19.6% (95%CrI: 2.6-54.4%) among MSM by 2030. CONCLUSIONS: New first-line antimicrobials for gonorrhoea treatment are needed. In the meantime, public health authorities should strengthen surveillance for AMR in N. gonorrhoeae and implement strategies for continued antimicrobial stewardship. Our models show the utility of long-term representative surveillance of gonococcal antimicrobial susceptibility data and can be adapted for use in, and for comparison with, other countries.

Early evidence of effectiveness of digital contact tracing for SARS-CoV-2 in Switzerland.

In the wake of the pandemic of coronavirus disease 2019 (COVID-19), contact tracing has become a key element of strategies to control the spread of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). Given the rapid and intense spread of SARS-CoV-2, digital contact tracing has emerged as a potential complementary tool to support containment and mitigation efforts. Early modelling studies highlighted the potential of digital contact tracing to break transmission chains, and Google and Apple subsequently developed the Exposure Notification (EN) framework, making it available to the vast majority of smartphones. A growing number of governments have launched or announced EN-based contact tracing apps, but their effectiveness remains unknown. Here, we report early findings of the digital contact tracing app deployment in Switzerland. We demonstrate proof-of-principle that digital contact tracing reaches exposed contacts, who then test positive for SARS-CoV-2. This indicates that digital contact tracing is an effective complementary tool for controlling the spread of SARS-CoV-2. Continued technical improvement and international compatibility can further increase the efficacy, particularly also across country borders.

Direct and indirect effects of the COVID-19 pandemic on mortality in Switzerland.

The direct and indirect impact of the COVID-19 pandemic on population-level mortality is of concern to public health but challenging to quantify. Using data for 2011-2019, we applied Bayesian models to predict the expected number of deaths in Switzerland and compared them with laboratory-confirmed COVID-19 deaths from February 2020 to April 2022 (study period). We estimated that COVID-19-related mortality was underestimated by a factor of 0.72 (95% credible interval [CrI]: 0.46-0.78). After accounting for COVID-19 deaths, the observed mortality was -4% (95% CrI: -8 to 0) lower than expected. The deficit in mortality was concentrated in age groups 40-59 (-12%, 95%CrI: -19 to -5) and 60-69 (-8%, 95%CrI: -15 to -2). Although COVID-19 control measures may have negative effects, after subtracting COVID-19 deaths, there were fewer deaths in Switzerland during the pandemic than expected, suggesting that any negative effects of control measures were offset by the positive effects. These results have important implications for the ongoing debate about the appropriateness of COVID-19 control measures.

Continuous population-level monitoring of SARS-CoV-2 seroprevalence in a large European metropolitan region.

Effective public health measures against SARS-CoV-2 require granular knowledge of population-level immune responses. We developed a Tripartite Automated Blood Immunoassay (TRABI) to assess the IgG response against three SARS-CoV-2 proteins. We used TRABI for continuous seromonitoring of hospital patients and blood donors (n = 72'250) in the canton of Zurich from December 2019 to December 2020 (pre-vaccine period). We found that antibodies waned with a half-life of 75 days, whereas the cumulative incidence rose from 2.3% in June 2020 to 12.2% in mid-December 2020. A follow-up health survey indicated that about 10% of patients infected with wildtype SARS-CoV-2 sustained some symptoms at least twelve months post COVID-19. Crucially, we found no evidence of a difference in long-term complications between those whose infection was symptomatic and those with asymptomatic acute infection. The cohort of asymptomatic SARS-CoV-2-infected subjects represents a resource for the study of chronic and possibly unexpected sequelae.

Impact of cross-border-associated cases on the SARS-CoV-2 epidemic in Switzerland during summer 2020 and 2021.

During the summers of 2020 and 2021, the number of confirmed cases of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections in Switzerland remained at relatively low levels, but grew steadily over time. It remains unclear to what extent epidemic growth during these periods was a result of the relaxation of local control measures or increased traveling and subsequent importation of cases. A better understanding of the role of cross-border-associated cases (imports) on the local epidemic dynamics will help to inform future surveillance strategies. We analyzed routine surveillance data of confirmed cases of SARS-CoV-2 in Switzerland from 1 June to 30 September 2020 and 2021. We used a stochastic branching process model that accounts for superspreading of SARS-CoV-2 to simulate epidemic trajectories in absence and in presence of imports during summer 2020 and 2021. The Swiss Federal Office of Public Health reported 22,919 and 145,840 confirmed cases of SARS-CoV-2 from 1 June to 30 September 2020 and 2021, respectively. Among cases with known place of exposure, 27% (3,276 of 12,088) and 25% (1,110 of 4,368) reported an exposure abroad in 2020 and 2021, respectively. Without considering the impact of imported cases, the steady growth of confirmed cases during summer periods would be consistent with a value of Re that is significantly above the critical threshold of 1. In contrast, we estimated Re at 0.84 (95% credible interval, CrI: 0.78-0.90) in 2020 and 0.82 (95% CrI: 0.74-0.90) in 2021 when imported cases were taken into account, indicating that the local Re was below the critical threshold of 1 during summer. In Switzerland, cross-border-associated SARS-CoV-2 cases had a considerable impact on the local transmission dynamics and can explain the steady growth of the epidemic during the summers of 2020 and 2021.

Historically High Excess Mortality During the COVID-19 Pandemic in Switzerland, Sweden, and Spain.

BACKGROUND: Excess mortality quantifies the overall mortality impact of a pandemic. Mortality data have been accessible for many countries in recent decades, but few continuous data have been available for longer periods. OBJECTIVE: To assess the historical dimension of the COVID-19 pandemic in 2020 for 3 countries with reliable death count data over an uninterrupted span of more than 100 years. DESIGN: Observational study. SETTING: Switzerland, Sweden, and Spain, which were militarily neutral and not involved in combat during either world war and have not been affected by significant changes in their territory since the end of the 19th century. PARTICIPANTS: Complete populations of these 3 countries. MEASUREMENTS: Continuous series of recorded deaths (from all causes) by month from the earliest available year (1877 for Switzerland, 1851 for Sweden, and 1908 for Spain) were jointly modeled with annual age group-specific death and total population counts using negative binomial and multinomial models, which accounted for temporal trends and seasonal variability of prepandemic years. The aim was to estimate the expected number of deaths in a pandemic year for a nonpandemic scenario and the difference in observed and expected deaths aggregated over the year. RESULTS: In 2020, the number of excess deaths recorded per 100 000 persons was 100 (95% credible interval [CrI], 60 to 135) for Switzerland, 75 (CrI, 40 to 105) for Sweden, and 155 (CrI, 110 to 195) for Spain. In 1918, excess mortality was 6 to 7 times higher. In all 3 countries, the peaks of monthly excess mortality in 2020 were greater than most monthly excess mortality since 1918, including many peaks due to seasonal influenza and heat waves during that period. LIMITATION: Historical vital statistics might be affected by minor completeness issues before the beginning of the 20th century. CONCLUSION: In 2020, the COVID-19 pandemic led to the second-largest infection-related mortality disaster in Switzerland, Sweden, and Spain since the beginning of the 20th century. PRIMARY FUNDING SOURCE: Foundation for Research in Science and the Humanities at the University of Zurich, Swiss National Science Foundation, and National Institute of Allergy and Infectious Diseases.

Assessing real-world vaccine effectiveness against severe forms of SARS-CoV-2 infection: an observational study from routine surveillance data in Switzerland.

BACKGROUND: In Switzerland, SARS-CoV-2 vaccination campaigns started in early 2021. Vaccine coverage reached 65% of the population in December 2021, mostly with mRNA vaccines from Moderna and Pfizer-BioNtech. Simultaneously, the proportion of vaccinated among COVID-19-related hospitalisations and deaths rose, creating some confusion in the general population. We aimed to assess vaccine effectiveness against severe forms of SARS-CoV-2 infection using routine surveillance data on the vaccination status of COVID-19-related hospitalisations and deaths, and data on vaccine coverage in Switzerland. METHODS: We considered all routine surveillance data on COVID-19-related hospitalisations and deaths received at the Swiss Federal Office of Public Health from 1 July to 1 December 2021. We estimated the relative risk of COVID-19-related hospitalisation or death for not fully vaccinated compared with fully vaccinated individuals, adjusted for the dynamics of vaccine coverage over time, by age and location. We stratified the analysis by age group and by calendar month. We assessed variations in the relative risk of hospitalisation associated with the time since vaccination. RESULTS: We included a total of 5948 COVID-19-related hospitalisations of which 1245 (21%) were fully vaccinated patients, and a total of 739 deaths of which 259 (35%) were fully vaccinated. We found that the relative risk of COVID-19 related hospitalisation was 12.5 (95% confidence interval [CI] 11.7-13.4) times higher for not fully vaccinated than for fully vaccinated individuals. This translates into a vaccine effectiveness against hospitalisation of 92.0% (95% CI 91.4-92.5%). Vaccine effectiveness against death was estimated to be 90.3% (95% CI 88.6-91.8%). Effectiveness appeared to be comparatively lower in age groups over 70 and during the months of October and November 2021. We also found evidence of a decrease in vaccine effectiveness against hospitalisation for individuals vaccinated for 25 weeks or more, but this decrease appeared only in age groups below 70. CONCLUSIONS: The observed proportions of vaccinated among COVD-19-related hospitalisations and deaths in Switzerland were compatible with a high effectiveness of mRNA vaccines from Moderna and Pfizer-BioNtech against hospitalisation and death in all age groups. Effectiveness appears comparatively lower in older age groups, suggesting the importance of booster vaccinations. We found inconclusive evidence that vaccine effectiveness wanes over time. Repeated analyses will be able to better assess waning and the effect of boosters.

Survival among people hospitalized with COVID-19 in Switzerland: a nationwide population-based analysis.

BACKGROUND: Increasing age, male sex, and pre-existing comorbidities are associated with lower survival from SARS-CoV-2 infection. The interplay between different comorbidities, age, and sex is not fully understood, and it remains unclear if survival decreases linearly with higher ICU occupancy or if there is a threshold beyond which survival falls. METHOD: This national population-based study included 22,648 people who tested positive for SARS-CoV-2 infection and were hospitalized in Switzerland between February 24, 2020, and March 01, 2021. Bayesian survival models were used to estimate survival after positive SARS-CoV-2 test among people hospitalized with COVID-19 by epidemic wave, age, sex, comorbidities, and ICU occupancy. Two-way interactions between age, sex, and comorbidities were included to assess the differential risk of death across strata. ICU occupancy was modeled using restricted cubic splines to allow for a non-linear association with survival. RESULTS: Of 22,648 people hospitalized with COVID-19, 4785 (21.1%) died. The survival was lower during the first epidemic wave than in the second (predicted survival at 40 days after positive test 76.1 versus 80.5%). During the second epidemic wave, occupancy among all available ICU beds in Switzerland varied between 51.7 and 78.8%. The estimated survival was stable at approximately 81.5% when ICU occupancy was below 70%, but worse when ICU occupancy exceeded this threshold (survival at 80% ICU occupancy: 78.2%; 95% credible interval [CrI] 76.1 to 80.1%). Periods with higher ICU occupancy (>70 vs 70%) were associated with an estimated number of 137 (95% CrI 27 to 242) excess deaths. Comorbid conditions reduced survival more in younger people than in older people. Among comorbid conditions, hypertension and obesity were not associated with poorer survival. Hypertension appeared to decrease survival in combination with cardiovascular disease. CONCLUSIONS: Survival after hospitalization with COVID-19 has improved over time, consistent with improved management of severe COVID-19. The decreased survival above 70% national ICU occupancy supports the need to introduce measures for prevention and control of SARS-CoV-2 transmission in the population well before ICUs are full.

Accounting for uncertainty during a pandemic.

We discuss several issues of statistical design, data collection, analysis, communication, and decision-making that have arisen in recent and ongoing coronavirus studies, focusing on tools for assessment and propagation of uncertainty. This paper does not purport to be a comprehensive survey of the research literature; rather, we use examples to illustrate statistical points that we think are important.

Regional excess mortality during the 2020 COVID-19 pandemic in five European countries.

The impact of the COVID-19 pandemic on excess mortality from all causes in 2020 varied across and within European countries. Using data for 2015-2019, we applied Bayesian spatio-temporal models to quantify the expected weekly deaths at the regional level had the pandemic not occurred in England, Greece, Italy, Spain, and Switzerland. With around 30%, Madrid, Castile-La Mancha, Castile-Leon (Spain) and Lombardia (Italy) were the regions with the highest excess mortality. In England, Greece and Switzerland, the regions most affected were Outer London and the West Midlands (England), Eastern, Western and Central Macedonia (Greece), and Ticino (Switzerland), with 15-20% excess mortality in 2020. Our study highlights the importance of the large transportation hubs for establishing community transmission in the first stages of the pandemic. Here, we show that acting promptly to limit transmission around these hubs is essential to prevent spread to other regions and countries.

Bayesian workflow for disease transmission modeling in Stan.

This tutorial shows how to build, fit, and criticize disease transmission models in Stan, and should be useful to researchers interested in modeling the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic and other infectious diseases in a Bayesian framework. Bayesian modeling provides a principled way to quantify uncertainty and incorporate both data and prior knowledge into the model estimates. Stan is an expressive probabilistic programming language that abstracts the inference and allows users to focus on the modeling. As a result, Stan code is readable and easily extensible, which makes the modeler's work more transparent. Furthermore, Stan's main inference engine, Hamiltonian Monte Carlo sampling, is amiable to diagnostics, which means the user can verify whether the obtained inference is reliable. In this tutorial, we demonstrate how to formulate, fit, and diagnose a compartmental transmission model in Stan, first with a simple susceptible-infected-recovered model, then with a more elaborate transmission model used during the SARS-CoV-2 pandemic. We also cover advanced topics which can further help practitioners fit sophisticated models; notably, how to use simulations to probe the model and priors, and computational techniques to scale-up models based on ordinary differential equations.

Socioeconomic position and the COVID-19 care cascade from testing to mortality in Switzerland: a population-based analysis.

BACKGROUND: The inverse care law states that disadvantaged populations need more health care than advantaged populations but receive less. Gaps in COVID-19-related health care and infection control are not well understood. We aimed to examine inequalities in health in the care cascade from testing for SARS-CoV-2 to COVID-19-related hospitalisation, intensive care unit (ICU) admission, and death in Switzerland, a wealthy country strongly affected by the pandemic. METHODS: We analysed surveillance data reported to the Swiss Federal Office of Public Health from March 1, 2020, to April 16, 2021, and 2018 population data. We geocoded residential addresses of notifications to identify the Swiss neighbourhood index of socioeconomic position (Swiss-SEP). The index describes 1·27 million small neighbourhoods of approximately 50 households each on the basis of rent per m2, education and occupation of household heads, and crowding. We used negative binomial regression models to calculate incidence rate ratios (IRRs) with 95% credible intervals (CrIs) of the association between ten groups of the Swiss-SEP index defined by deciles (1=lowest, 10=highest) and outcomes. Models were adjusted for sex, age, canton, and wave of the epidemic (before or after June 8, 2020). We used three different denominators: the general population, the number of tests, and the number of positive tests. FINDINGS: Analyses were based on 4 129 636 tests, 609 782 positive tests, 26 143 hospitalisations, 2432 ICU admissions, 9383 deaths, and 8 221 406 residents. Comparing the highest with the lowest Swiss-SEP group and using the general population as the denominator, more tests were done among people living in neighbourhoods of highest SEP compared with lowest SEP (adjusted IRR 1·18 [95% CrI 1·02-1·36]). Among tested people, test positivity was lower (0·75 [0·69-0·81]) in neighbourhoods of highest SEP than of lowest SEP. Among people testing positive, the adjusted IRR was 0·68 (0·62-0·74) for hospitalisation, was 0·54 (0·43-0·70) for ICU admission, and 0·86 (0·76-0·99) for death. The associations between neighbourhood SEP and outcomes were stronger in younger age groups and we found heterogeneity between areas. INTERPRETATION: The inverse care law and socioeconomic inequalities were evident in Switzerland during the COVID-19 epidemic. People living in neighbourhoods of low SEP were less likely to be tested but more likely to test positive, be admitted to hospital, or die, compared with those in areas of high SEP. It is essential to continue to monitor testing for SARS-CoV-2, access and uptake of COVID-19 vaccination and outcomes of COVID-19. Governments and health-care systems should address this pandemic of inequality by taking measures to reduce health inequalities in response to the SARS-CoV-2 pandemic. FUNDING: Swiss Federal Office of Public Health, Swiss National Science Foundation, EU Horizon 2020, Branco Weiss Foundation.

Estimating the cumulative incidence of SARS-CoV-2 with imperfect serological tests: Exploiting cutoff-free approaches.

Large-scale serological testing in the population is essential to determine the true extent of the current SARS-CoV-2 pandemic. Serological tests measure antibody responses against pathogens and use predefined cutoff levels that dichotomize the quantitative test measures into sero-positives and negatives and use this as a proxy for past infection. With the imperfect assays that are currently available to test for past SARS-CoV-2 infection, the fraction of seropositive individuals in serosurveys is a biased estimator of the cumulative incidence and is usually corrected to account for the sensitivity and specificity. Here we use an inference method-referred to as mixture-model approach-for the estimation of the cumulative incidence that does not require to define cutoffs by integrating the quantitative test measures directly into the statistical inference procedure. We confirm that the mixture model outperforms the methods based on cutoffs, leading to less bias and error in estimates of the cumulative incidence. We illustrate how the mixture model can be used to optimize the design of serosurveys with imperfect serological tests. We also provide guidance on the number of control and case sera that are required to quantify the test's ambiguity sufficiently to enable the reliable estimation of the cumulative incidence. Lastly, we show how this approach can be used to estimate the cumulative incidence of classes of infections with an unknown distribution of quantitative test measures. This is a very promising application of the mixture-model approach that could identify the elusive fraction of asymptomatic SARS-CoV-2 infections. An R-package implementing the inference methods used in this paper is provided. Our study advocates using serological tests without cutoffs, especially if they are used to determine parameters characterizing populations rather than individuals. This approach circumvents some of the shortcomings of cutoff-based methods at exactly the low cumulative incidence levels and test accuracies that we are currently facing in SARS-CoV-2 serosurveys.

Estimation of SARS-CoV-2 mortality during the early stages of an epidemic: A modeling study in Hubei, China, and six regions in Europe.

BACKGROUND: As of 16 May 2020, more than 4.5 million cases and more than 300,000 deaths from disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have been reported. Reliable estimates of mortality from SARS-CoV-2 infection are essential for understanding clinical prognosis, planning healthcare capacity, and epidemic forecasting. The case-fatality ratio (CFR), calculated from total numbers of reported cases and reported deaths, is the most commonly reported metric, but it can be a misleading measure of overall mortality. The objectives of this study were to (1) simulate the transmission dynamics of SARS-CoV-2 using publicly available surveillance data and (2) infer estimates of SARS-CoV-2 mortality adjusted for biases and examine the CFR, the symptomatic case-fatality ratio (sCFR), and the infection-fatality ratio (IFR) in different geographic locations. METHOD AND FINDINGS: We developed an age-stratified susceptible-exposed-infected-removed (SEIR) compartmental model describing the dynamics of transmission and mortality during the SARS-CoV-2 epidemic. Our model accounts for two biases: preferential ascertainment of severe cases and right-censoring of mortality. We fitted the transmission model to surveillance data from Hubei Province, China, and applied the same model to six regions in Europe: Austria, Bavaria (Germany), Baden-Württemberg (Germany), Lombardy (Italy), Spain, and Switzerland. In Hubei, the baseline estimates were as follows: CFR 2.4% (95% credible interval [CrI] 2.1%-2.8%), sCFR 3.7% (3.2%-4.2%), and IFR 2.9% (2.4%-3.5%). Estimated measures of mortality changed over time. Across the six locations in Europe, estimates of CFR varied widely. Estimates of sCFR and IFR, adjusted for bias, were more similar to each other but still showed some degree of heterogeneity. Estimates of IFR ranged from 0.5% (95% CrI 0.4%-0.6%) in Switzerland to 1.4% (1.1%-1.6%) in Lombardy, Italy. In all locations, mortality increased with age. Among individuals 80 years or older, estimates of the IFR suggest that the proportion of all those infected with SARS-CoV-2 who will die ranges from 20% (95% CrI 16%-26%) in Switzerland to 34% (95% CrI 28%-40%) in Spain. A limitation of the model is that count data by date of onset are required, and these are not available in all countries. CONCLUSIONS: We propose a comprehensive solution to the estimation of SARS-Cov-2 mortality from surveillance data during outbreaks. The CFR is not a good predictor of overall mortality from SARS-CoV-2 and should not be used for evaluation of policy or comparison across settings. Geographic differences in IFR suggest that a single IFR should not be applied to all settings to estimate the total size of the SARS-CoV-2 epidemic in different countries. The sCFR and IFR, adjusted for right-censoring and preferential ascertainment of severe cases, are measures that can be used to improve and monitor clinical and public health strategies to reduce the deaths from SARS-CoV-2 infection.

Pattern of early human-to-human transmission of Wuhan 2019 novel coronavirus (2019-nCoV), December 2019 to January 2020.

Since December 2019, China has been experiencing a large outbreak of a novel coronavirus (2019-nCoV) which can cause respiratory disease and severe pneumonia. We estimated the basic reproduction number R0 of 2019-nCoV to be around 2.2 (90% high density interval: 1.4-3.8), indicating the potential for sustained human-to-human transmission. Transmission characteristics appear to be of similar magnitude to severe acute respiratory syndrome-related coronavirus (SARS-CoV) and pandemic influenza, indicating a risk of global spread.